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Purpose: Manipulation of the radio-sensitivity of the nucleotide-base driven by the spin blockade 
mechanism of diffusive free radicals against ionizing radiation. 

Materials and methods: We theoretically propose a mechanism which uses the simultaneous 
application of circularly polarized light and an external magnetic field to control the polarization 
of the free radicals and create S=l electron-hole spin excitations (excitons) on nucleotide-base. We 
deploy an ab-initio molecular dynamics model to calculate the characteristic parameters of the light 
needed for optical transitions. 

Results: As a specific example, we present the numerical results calculated for a Guanine, in the 
presence of an OH free radical. To increase the radio-resistivity of this system, a blue light source 
for the optical pumping and induction of excitons on guanine can be used. 

Conclusions: The effect of spin-injection on the formation of a free energy barrier in diffusion 
controlled chemical reaction pathways leads to the control of radiation-induced base damage. The 
proposed method allows us to manipulate and partially suppress the damage induced by ionizing 
radiation. 



PACS numbers: 

INTRODUCTION 

Ionizing radiation is both hazardous and beneficial to 
living organisms, and is extensively used for cancer treat- 
ment in radiation therapy [Ij . A major problem in the ap- 
plication of ionizing radiation to cancer treatment is the 
protection of normal cells and tissues against unavoid- 
able exposure to radiation during radiation treatment. 
It is now well understood that the ionization or excita- 
tion of the DNA molecules, either directly or indirectly, 
can lead to DNA single or double strand breaks. As a 
result, misrepaired DNA molecules can lead to specific 
genetic aberrations and/or mutations which could cause 
carcinogenesis in normal cells or lead to fatal damage in 
normal or cancer cells E]. It has been shown that 
low linear-energy-transfer (LET) ionizing radiation cre- 
ates approximately 1,000 single strand breaks (SSBs) and 
40 double strand breaks (DSBs) per gray (lGy=lJ/Kg) 
in typical mammalian cells |H [5) [7] • The level of DNA 
molecular base damage is around 2,500 to 25,000 per Gy 
in a cell, which is about 2.5 to 25 times the yield of sugar- 
phosphate induced damage in the DNA backbone [U E] • 
In indirect mechanisms, the water molecules surround- 
ing the DNA molecule which compose 80% of a cell, may 
be excited by ionizing radiation in form of free radicals, 
e.g., a charged neutral hydroxyl (OH). The motion of 
OH-radicals which are randomly produced throughout 
the cell is governed by diffusion processes. Massive DNA 
damage can result from a large number of DNA dehy- 
drogenations caused by free radicals. For example, a free 
radical can diffuse to reach a DNA molecule and remove 
a hydrogen ion from it to form a water molecule. De- 
tailed studies at the molecular level is necessary to bring 
the radiation-induced DNA damage under control. 



METHOD 

In this work, we apply a quantum physical description 
of molecular interactions to propose a mechanism that 
could allow the manipulation of DNA radio-sensitivity. 
In particular the Pauli exclusion principle [5] which pre- 
vents two electrons with parallel spin form occupying a 
single spatial orbital, plays a major role and is used to 
magnetically manipulate the diffusion of hydroxyl radi- 
cals and the OH-DNA relative motion. It has been shown 
in studies in semiconductor physics and quantum optics 
that the Pauli exclusion principle can be used to rec- 
tify electrical currents passing through weakly coupled 
quantum dots |TTJ and to induce ferromagnetic ordering 
by photo-generated carriers in magnetic semiconductor 
hetero-structures [T2"] . 

A free radical carries an odd number of electrons with 
an unpaired spin in the outermost open shell. Due to 
the reduction of the exchange interaction, the pairing of 
opposite spin electrons in the open orbital of the free 
radical with an electron in a DNA molecule makes free 
radicals highly reactive. In the process of dehydrogena- 
tion of a DNA molecule by free radicals an unpaired 
hole (a half-empty orbital) is transferred to the DNA. 
In the absence of spin-orbit coupling and hyperfine inter- 
action the spin of transfered electron is conserved. The 
electronic ground state of DNA-molecule is S — spin- 
singlet (in the absence of an external magnetic field) . The 
OH-radical which contains nine electrons is a doubly de- 
generate ground state with S z — ±1/2, where we have 
conveniently taken the quantization axis along the z-axis. 
The degeneracy of the ground state can be lifted by ap- 
plying a weak magnetic field which couples to the electron 
spin through the Zeeman interaction [S], Ez = ffMsS 1 - B. 
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Here Ez is the Zeeman energy, g is the electron <?-factor 
(g w 2), /xg is the Bohr magneton (/ie w 5.8 x 10~ 5 
eV/Tesla), and _B is the strength of external magnetic 
field. 

In a random interaction of radiation with a biological 
system the initial direction of OH-radical magnetic mo- 
ment immediately after its generation is also random. 
However, by applying a weak external magnetic field 
(-Bext) (which defines the quantization axis) and using 
a circularly polarized light field parallel to the direction 
of the light propagation, as shown in Fig jlji, a molecular 
transition corresponding to A J = ±1 can be induced by 
means of optical pumping (T3] of the OH-radicals [H] . 
Here J denotes the total angular momentum of diatomic 
OH-radical [15 . Alternatively, tchniques such as electron 
spin resonance (ESR) can be used to achieve strong po- 
larization of free radicals, as recent advances in ESR have 
demonstrated the capability of detecting the transfer of 
electron spin polarization between radicals [TU] . In this 
case microwaves can be used for the optical transitions. 
In a similar fashion, by applying a second circularly po- 
larized light field one may excite an electron-hole pair 
(exciton) in the DNA molecule. Because the circularly 
polarized light carries angular momentum ±1, the exci- 
ton has a particular spin polarization. Here the spin of 
exciton is S = 1 with polarization along the light propa- 
gation direction (because of angular momentum selection 
rules) . Fig. [I] schematically shows the generation of the 
optically pumped exciton by circularly polarized light. 
The injection of photo-electrons with the spin out of equi- 
librium may lead to a dramatic effect in the collective dy- 
namical behavior of DNA-molecules and the interaction 
with OH-radicals. For example the OH-DNA repulsive 
magnetic force provides a potential barrier which blocks 
the diffusion pathway (see Fig. [lj>) of OH-radicals to- 
ward the DNA-molecules. This is expected to hinder the 
DNA dehydrogenation and consequently increase the cell 
radio-resistivity. To verify this hypothesis, an a&-initio 
molecular dynamical model, which is the mathematical 
formulation that governs the appropriate dynamics of the 
molecular system [IB] is deployed. We have used the Car- 
Parrincllo molecular dynamics (CPMD) [18, 19J model, 
in which the potential energy of the system can be calcu- 
lated on-the-fly, as needed for the conformations of the 
dynamical trajectory to simulate the chemical reaction 
pathways. Because the absorption of a circularly polar- 
ized photon alters the local electronic state of a DNA- 
molecule, we confine our simulation to a particular seg- 
ment, e.g., only a part of the DNA where the injected 
exciton is localized and the optical transition takes place. 
To illustrate this, let us consider a system of interest con- 
sisting of a DNA nucleotide base, (e.g., guanine) in the 
presence of the OH-radical. We assume that a photon 
with circular polarization interacting with guanine can 
induce an optical transition in the form of an S = 1 
exciton. Here we investigate the effect of an exciton pro- 



duced in this way on the guanine-dehydrogenation path- 
way, assuming that another photon generated through 
interactions with ionizing radiation (such as radiotherapy 
x-rays or cosmic rays) creates a free radical in the vicin- 
ity of guanine. Because the local density of free radicals 
and excitons are large and are comparable, the events 
described in our calculation can be observed with rea- 
sonable probability. We adopt computational parameters 
and variables needed for the CPMD calculation of the dy- 
namical trajectory of the gas phase nucleotide bases in 
the presence of OH-radicals following Refs. [17 , where 
the consistency of CPMD results for guanine with other 
quantum chemistry approaches has been investigated. 

RESULTS 

We identify the dehydrogenation of the nucleotide 
bases as a function of their spin multiplicity. The ground 
and excited states of the nucleotide correspond to spin 
singlet (S = 0), and spin triplet (S — 1) states. The 
latter can be realized through the application of circu- 
larly polarized light as discussed above (see Fig. [TJ. 
Our CPMD is implemented in a plane-wave basis within 
local spin density approximation (LSDA) with an en- 
ergy cutoff of 70 Rydberg (Ry), and with Becke [20] 
exchange and Lee- Yang-Parr (BLYP) gradient-corrected 
functional [3T]. Norm conserving ultrasoft Vanderbilt 
pseudo-potentials were used for oxygen, hydrogen, nitro- 
gen and carbon. The CPMD micro-canonical dynamics 
(constant energy ensemble) were performed after wave- 
function optimization following dynamical equilibration 
at T=300K and re-quenching of the wave-function. An 
isolated cubic cell of length 13.229 A with Poisson solver 
of Martyna and Tuckerman (22] was used. Our CPMD 
studies consist of two classes of spin-restricted calcula- 
tions, as the total spin along the quantum axis is sub- 
jected to the constraints S z = 1/2, and 3/2, correspond- 
ing to doublet and quartet spin configurations. In both 
calculations the initial distance between OH-radical and 
nucleotide is considered to be about 1.5 A. We selectively 
choose an initial coordinate for OH-radical in the neigh- 
borhood of the nucleotide where the Hydrogen transfer 
shows a reactive path in normal state of DNA (the dou- 
blet spin configuration in the absence of circularly polar- 
ized light and magnetic field). 

The initial and final states of the molecules are shown 
in Figs. |2]4j The final configurations of the molecules 
have been obtained after 0.6 ps where the rearrangement 
of the atomic coordinates have been deduced from a dy- 
namical trajectory calculated by CPMD. According to 
our results, a rapid dehydrogenation of the nucleotides 
takes place for a system with S z = 1/2 (total spin- 
doublet) as shown in Fig. [3] This process leads to the 
formation of a water molecule. In contrast, as shown 
in Fig. [4] in the quartet spin configuration the repul- 
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FIG. 1: Schematically shown in (a) the injection of photo- 
generated electrons in DNA-nucleotides with spin polarization 
(shown by arrows) along the direction of circularly polarized 
light and external magnetic field. The net magnetic force be- 
tween two parallel magnetic moments localized in OH and 
DNA-nucleotide is repulsive. This is similar to two separated 
magnetic moments which interact like Heisenberg antiferro- 
magnetic exchange coupling (b). 



FIG. 3: The state of de-hydrogenated Guanine by OH free 
radical at t = 0.6 ps. The polarization state of the system is 
spin doublet (S = 1/2). 




FIG. 2: Initial state of Guanine molecule in the presence of 
irradiated induced OH free radical. 



sive exchange interaction, analogous to Heisenberg anti- 
ferromagnetic coupling which originates from the Pauli 
exclusion principle, blocks the exchange of hydrogen and 
hence the chemical reaction. In Fig. [5] the evolution of 
the Ni Hydrogen in the guanine and free radical oxygen 
distance is shown. As it is seen the abstraction of Hy- 
drogen occurs around t sa 50fs in the spin singlet state of 
guanine, and the injection of S — 1 exciton in guanine 
blocks the hydrogen abstraction. Fig. [6]shows the Kohn- 
Sham energies (equivalent to potential energy in classical 
molecular dynamics) of the spin singlet and spin triplet 
of the Guanine in the presence of the OH free radical as 
a function of time, calculated by the CPMD at T=300K 
corresponding to a canonical dynamics (constant temper- 
ature ensemble). A drop in Kohn-Sham energy in spin 
singlet multiplicity is indication of dehydrogenation of 
Hjvi in the guanine by OH free radical. To systematically 
check the convergence of the results, we increased the 
size of the molecule by adding sugar-phosphate rings to 
guanine and found that this has no influence on the spin- 
blocking effect. To estimate the energy needed for the po- 
larization of the nucleotide in the absence of OH-radicals, 
we calculated the energy of the ground and excited states 
of the gas-phase nucleotide in spin singlet and triplet mul- 
tiplicities. For guanine we calculated the spin singlet-spin 
triplet energy gap A = -^triplet - -Esingiet ~ 2.68 eV. This 
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FIG. 5: The evolution of the distances in Angstrom from oxy- 
gen atom in the OH radical to the Hjvi in the guanine as a 
function of guanine spin multiplicity. The hydrogen abstrac- 
tion occurs around t ~ 50fs in spin singlet state of guanine. 
The injection of S = 1 exciton in guanine blocks the hydrogen 
abstraction. 



FIG. 4: The state of radio-resistive Guanine at t — 0.6 ps. 
The polarization state of the system is spin quartet (S = 3/2) 
induced by circularly polarized light in the presence of weak 
magnetic field. Due to injected polarized photo-electrons lo- 
calized in Guanine, the dehydrogenated Guanine does not 
form. 



provides an estimate for the frequency of the circularly 
polarized light, which is within the range of the visible 
spectrum of the electromagnetic waves, A = 463 nm (light 
blue). To calculate the stored magnetic energy due to 
the optical injection of spin, we calculated the energy of 
the gas-phase nucleotide in the presence of one OH-frcc 
radical with spin doublet and quartet multiplicities. For 
the molecules shown in Fig. [2j we find the energy gap 



Ai = E, 



quartet 



E, 



doublet 



3.54 eV. Here the excessive 



magnetic energy which originated from spin-spin repul- 
sive interactions (which resemble the anti-ferromagnetic 
exchange interaction in the Heisenberg model) can be 
deduced to be Ai — Ao ~ 0.86 eV. This energy can be 
interpreted as the excessive energy barrier due to the 
alignment of the spins in the DNA molecule and OH, 
and is the source of the magnetic repulsive force which 
makes the diffusion of OH toward DNA-molecules less 
likely. This is in agreement with the results obtained 
from CPMD, shown in Figs. 2]|4 In addition, by switch- 
ing the polarization of one of the light sources to the 
opposite direction, the relative direction of the DNA-OH 
polarization switches to antiparrallel, and hence the mag- 
netic repulsive force changes to an attractive force that 
lowers the OH diffusion barrier and decreases the radio- 
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FIG. 6: The Kohn-Sham energy as a function of time and 
spin multiplicity of Guanine, spin triplet (top) and spin sin- 
glet (bottom). A drop in Kohn-Sham energy in spin singlet 
multiplicity is indication of dehydrogenation of Hjvi in the 
guanine by OH free radical. 



sensitivity of the DNA-molecule. 

After photon absorption, the nucleotide is spin polar- 
ized along the direction determined by the polarization 
state and the propagation direction of the circularly po- 
larized light. The polarized state of the nucleotide then 
decays quantum mechanically to its unpolarized ground 
state either by spontaneous photon emission (electron- 
hole recombination) or through photon-electron spin de- 
coherence. There are radiative and non-radiative chan- 
nels that contribute to this process. Since spin-orbit 
coupling governs one of the decay mechanisms in non- 
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radiative channel, we use Fermi's golden rule to estimate 
the life-time of the triplet state. It then follows that 
T T ^ S = (2ir/h 2 )VL J d 3 k/(2ir) 3 \(T\Hso\S)\ 2 6(uJo — w) 
[5]. Here, T^^s is the transition rate from the spin- 
triplet (T) to the spin-singlet (S) state, f2 is the vol- 
ume, k — 27r/A is the emitted photon wave-number, 
Hso = -e/(2ra 2 c 2 )EiI 1 s , i • (ft x V^ K s(n)) is the 
spin-orbit Hamiltonian, m is the electron mass, c is the 
speed of light, Si, pi are the spin and momentum of 
the ith electron, $ks is the Kohn-Sham effective po- 
tential, and uoq = Aq /h. This calculation shows that 



techniques for radiation therapy and radiation protection 
purposes. 
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100 ps. However, the electronic relaxed 
excitonic states with emperical lifetime of several 100 ps 
has been reported recently [23J . The spin-triplet life time 
of nucleotide r turns out to be significantly larger than 
the dehydrogenation time scale. It is therefore possible to 
increase the radio-resistivity of the DNA molecule within 
this time scale through optically pumped spin polariza- 
tion. It is important to compare r with other time-scales 
in the process. The initial ionization takes place in about 
1 fs (10 -15 second). The primary free radicals produced 
by ejection of electrons have a life time of nearly 100 ps, 
and the reported OH radical life-time is about 1 ns (10 -9 
second) |3J. The electron spin-lattice relaxation time of 
the OH radical has been estimated to be approximately 
between 0.1 and 0.5 ns in water at room temperature |10j . 

In order to estimate the technical requirements for the 
above described approach one could assume that an aque- 
ous solution of DNA will be irradiated with a dose of lGy 
(U/Kg). It is known [24] that lOOeV of absorbed pho- 
ton/electron energy produce about 6 OH radicals. There- 
fore lGy of radiation produce about 4 x 10 13 OH radicals 
in 0.1 cm 3 of water. If the number of injected excitons can 
be exceeded up to at least ten times, by applying a laser 
pump with moderate intensity it is possible to increase 
significantly the resistance of DNA-molecules against ir- 
radiation. For example at a dose rate of 1.4 Gy/min a 
laser pump power of P = 10NoH^o/ T ir ~ 12 x 10~ 6 
watt would be required, which is well within the techni- 
cally achievable limits. 



CONCLUSION 

In conclusion, we have theoretically explored a mech- 
anism which involves the injection of spin polarized ex- 
citons in DNA molecules to control and manipulate the 
radio-sensitivity of cells by using a circularly polarized 
light field and external magnetic field. The mechanism 
proposed here is based on the selection rules applica- 
ble to optical transitions between energy levels of the 
DNA-molecules and optical pumping of the OH-radicals, 
and we have employed a microscopic ab-initio molecular 
dynamics model to computationally study the dehydro- 
genation mechanism at the molecular level. The results 
of this study may be used as a guideline to develop new 



[9 

[io; 
in 

[12 
[13 

[14; 

[15 
[16 

[17 

[is; 

[19 



Faiz M. Khan, The Physics of Radiation Therapy, 

(Williams & Wilkins, Baltimore, 1994). 

L. G. Burdelya, V. I. Krivokrysenko, T. C. Tallant, E. 

Strom, A. S. Gleiberman, D. Gupta, O. V. Kurnasov, F. 

L. Fort, A. L. Osterman, J. A. DiDonato, E. Feinstein, 

and A. V. Gudkov, Science 320 226-230 (2008). 

Eric J. Hall, Radiobiology for the Radiologist, (Lippincott 

Williams & Wilkins, Baltimore, Fifth Edition, 2000). 

J.F. Ward, Prog. Nucleic Acid Res. Mol. Biol., 35 95 

(1988); Radiat. Res. 142, 362 (1995); Errata, Radiat. 

Res. 143, 355 (1995). 

D. T. Goodhead, Int. J. Radiat. Biol. 65, 7 (1994). 
H. Nikjoo et ai, Int. J. Radiat. Bio. 71, 467 (1997). 
V. A. Semenenko, and R. D. Stewart, Radiat. Res. 161, 
451 (2004). 

L. D. Landau, and E. M. Lifshitz, Quantum Mechanics: 
Non-Relativistic Theory (Pergamon, Oxford, 2003). 
D. M. Bartels, A. D. Trifunac, R. G. Lawler, Chem. Phys. 
Lett. 152, 109 (1988); William S. Jenks and Nicholas J. 
Turro, J. Am. Chem. Soc. 112 9009 (1990). 
N. C. Verma and R. W. Fessenden, J. Chem. Phys. 65, 
2139 (1976); B. Brocklehurst, J. Chem. Soc. Faraday 
Trans. II 75, 123 (1979); B. Bhattacharjee, and Ranjan 
Das, Molecular Physics, 105, 1053 (2007). 
K. Ono, K., D. G. Austing, Y. Tokura, S. Tarucha, Sci- 
ence 297 1313 (2002). 

S. Koshihara, A. Oiwa, M. Hirasawa, S. Katsumoto, Y. 
Iye, C. Urano, H. Takagi, and H. Munekata, Phys. Rev. 
Lett. 78, 4617 (1997); I. Zutic, J. Fabian, and S. Das 
Sarma, Rev. Mod. Phys. 76, 323 (2004). 
Thad G. Walker and William Happer, Rev. Mod. Phys. 
69, 629 (1997); William Happer, Rev. Mod. Phys. bf 44, 
169 (1972). 

G.H. Dieke, and H. M. Crosswhite, J. Quant. Spectrosc. 
Radiat. Transf. 2, 97 (1962); Sebastiaan Y. van de Meer- 
akker, Nicolas Vanhaecke, Mark P. van der Loo, Gerrit 

C. Groenenboom, and Gerard Meijer, Phys. Rev. Lett. 
95, 013003 (2005). 

G. Herzberg, The Spectra and Structures of Simple Free 
Radicals (Dover Publications Inc. New York, 1971). 

D. Marx, D., and J. Hutter, Ab initio molecular dynam- 
ics: Theory and Implementation, Modern Methods and 
Algorithms of Quantum Chemistry, J. Grotendorst (Ed.), 
(John von Neumann Institute for Computing, Jiilich, 
NIC Series, Vol. 1, page 301-449, 2000). 

C. J. Mundy, M. E. Colvin, A. A. Quong, J. Phys. Chem. 

A. 106 10063 (2002); Y. Wu, C. J. Mundy, M. E. Colvin, 

R. Car, J. Phys. Chem. A. 108 2922 (2004). 

R. Car, and M. Parrinello, Phys. Rev. Lett. 55, 2471 

(1985). 

J. Hutter, P. Ballone, M. Bernasconi, P. Focher, E. Fois, 
S. Goedecker, M. Parrinello, M. E. Tuckerman, CPMD 
code, version 3.13, MPI fuer Festkoerperforschung, 
Stuttgart IBM Zurich Research Laboratory, 1990-2008. 



[20] A. D. Becke, Phys. Rev. A 38, 3098 (1988). 
[21] C. Leo, W. Yang, and R. G. Parr, Phys. Rev. B 37, 785 
(1988). 

[22] G. J. Martyna, M. E. Tuckerman, J. Chem. Phys. 110, 
2810 (1999). 

[23] I. Buchvarov, Q. Wang, M. Raytchev, A. Trifonov, T. 

Fiebig, Proc. Natl. Acad. Sci. 104, 4794 (2007); Q. 

Wang, and T. Fiebig in Charge Migration in DNA, Ed. 

T. Chakraborty (Springer- Verlag, Berlin, 2007). 
[24] Yamaguchi Hiroshi et al, J. Radiat. Res. 46, 333 (2005). 



